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ABSTRACT 

We discuss the reference time to of afterglow light curves in the context of 
the standard internal-external shock model. The decay index of early afterglow 
is very sensitive to the reference time one chooses. In order to understand the 
nature of early afterglow, it is essential to take a correct reference time. Our 
simple analytic model provides a framework to understand special relativistic ef- 
fects involved in early afterglow phase. We evaluate light curves of reverse shock 
emission as well as those of forward shock emission, based on full hydrodynamic 
calculations. We show that the reference time does not shift significantly even in 
the thick shell case. For external shock emission components, measuring times 
from the beginning of the prompt emission is a good approximation and it does 
not cause an early steep decay. In the thin shell case, the energy transfer time 
from fireball ejecta to ambient medium typically extends to thousands of sec- 
onds. This might be related to the shallow decay phases observed in early X-ray 
afterglow at least for some bursts. 

Subject headings: gamma rays: bursts — relativity — hydrodynamics 



1. Introduction 

It is well known that gamma-ray burst (GRB) afterglow decays as a power-law L oc 
(t — to)~ a . The temporal decay index a, together with the spectral index, provides us 
precious information about GRB jets and their environment. In the pre-Swift era afterglow 
observations start typically a few hours after a burst. In such a late phase, the decay index 
is insensitive to the choice of the reference time to, and GRB trigger time is often used in 
afterglow modelings. 

The multi-wavelength observatory Swift was launched in Nov 2004. Thanks to its fast 
pointing capabilities Swift is disclosing the early afterglow phase. One of unexpected finds 
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by Swift is that early X-ray afterglows show a canonical behavior, where light curves include 
three components: (1) a steep decay component, (2) a shallow decay component and (3) a 
"normal" decay component. On top of this canonical behavior, many events have superim- 
posed X-ray flares (Zhang et al. 2006a; Nousek et al. 2006; Chincarini et al 2005; O'Brien 
et al. 2006). The transition from the early steep decay to the shallow decay typically oc- 
curs at several hundred seconds, and the timescale is comparable to the duration of rather 
long GRBs. When discussing the early afterglow and its connection to the prompt emission 
component, the decay index is very sensitive to the reference time t one chooses. Correctly 
choosing t is therefore essential to derive the right index as well as to interpret each com- 
ponent in the canonical light curve (Piro et al. 2005; Tagliaferri et al. 2005; Quimby et al. 
2006). 

Tagliaferri et al. (2005) investigated the first two bursts GRB 050126 and GRB 050219a 
which have an X-ray light curve well sampled by the X-Ray Telescope on board Swift. They 
sought for a possible delay of the afterglow onset by fitting the early X-ray light curves (the 
components 1 and 2 which we have discussed above) with a single power-law model. In 
both cases, the decaying light curves can be fitted if the onset of the afterglow is shifted to 
to ~ 100 s after the burst trigger with a single power-law. However, while in the case of 
GRB 050126 the light curve does not allow us to clearly state whether a broken power-law 
modeling is better than a single power-law model, for GRB 050219a a broken power- law 
definitively provides a better fit. 

In the standard GRB model the time shift between the GRB trigger and the reference 
time to is expected to be "small" . The early steep decay should not be an artifact due to a 
wrong choice of to- Lazzati and Begelman (2006) studied forward shock emission, based on a 
simple energy injection model. Their numerical light curves show that measuring times from 
the beginning of the prompt phase is a good approximation. The early steepening and X-ray 
flares are likely to be produced by another mechanism (e.g. internal shocks; Burrows et al. 
2005; Falcone et al. 2006; Zhang et al. 2006a; Nousek et al. 2006; Ioka, Kobayashi & Zhang 
2005; Fan & Wei 2005). Recently long-lasting soft emission is reported in a short burst 
GRB 050724 (Barthelmy et al. 2005). Such a soft component was hinted at in the sum of 
multiple short BATSTE GRBs, and it might be the onset of short burst afterglow (Lazzati, 
Ramirez- Ruiz & Ghisellini 2001). It is therefore of interest to quantitatively examine the t 
issue. 

In this paper, we study the physics and timescales involved in early afterglow stage, 
and give more direct and clear arguments for the afterglow slopes. In §2 we study a simple 
analytic model. In §3 we evaluate light curves of reverse shock emission as well as that of 
forward shock emission, based on full hydrodynamic calculations. In §4 we address how 
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inhomogeneity of a fireball affects an early afterglow light curve. Conclusions and discussion 
are given in §5. 



Let R be a radius of a forward shock expanding with a Lorentz factor r » 1 into 
homogeneous ambient medium. Since the shock is moving toward us at almost the speed of 
light c, the difference of the the observed times between a photon emitted at R and another 
emitted at R + dR is dt ~ dR/2cT 2 . Although the origin of observed time is arbitrary, a 
natural definition of observed time is given by the delay of photons emitted from a shock 
front at a lab time t with respect to the photons emitted from the "explosion" at R = and 
i = 0. The dashed line in figure 1 depicts the trajectory of the photon from the explosion. 
In this paper, i and t denote the lab and observed time since the explosion, respectively. 
Distance, velocity and the corresponding Lorentz factors are measured in the lab frame. 
Thermodynamic quantities (pressure and density) are measured in the local fluid frame. 
The observed time is given by 



The shock radius is almost proportional to the lab time, i ~ R/c. Considering that internal 
shocks occur and produce gamma-rays at radii much smaller than the deceleration radius of 
the fireball, the prompt gamma-rays associated with the outermost element of the fireball 
should propagate practically on the dashed line in figure 1. The GRB trigger almost coincides 
with the explosion t = 0. The difference is small, and order of the variability timescale of 
the prompt emission, because the variability timescale directly reflects the inhomogeneity 
scale in a fireball (Kobayashi, Piran & Sari 1997). The dashed line in figure 1 will be called 
the gamma-ray front in the following sections. 

A fireball with an initial Lorentz factor To decelerates when it collects a large volume of 
ambient material with a mass density p\. Equalizing energy of the shocked ambient material 
47ri? 3 pic 2 ro/3 and fireball energy E, we obtain the deceleration radius Rd = l/F^ 3 where 
I = (3i?/47rp 1 c 2 ) 1 / 3 is the Sedov length. Since the fireball density decreases as it expands, 
there is a possibility that a reverse shock evolves from Newtonian to relativistic during the 
propagation. A relativistic reverse shock reduces considerably the Lorentz factor of the 
fireball material which it crosses. In such a case, the energy of the shocked ambient medium 
is still negligible at R = l/T^ 2 , because we have assumed r ~ r at the deceleration to 
estimate Rd- The reverse shock crosses the fireball shell at R = Z 3 / 4 A 1 / 4 and all the fireball 
material decelerates where A is a fireball shell width (see Sari & Piran 1995 and Kobayashi, 
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Piran & Sari 1999 for the details). In summary, the evolution of fireballs are classified into 
two cases depending on the value of A relative to a critical value A c = l/T^/ 3 (Sari & 
Piran 1995). If A is smaller than the critical value, it is called the thin shell case 1 . The 
reverse shock is always in the Newtonian regime, and it is too weak to slow down the fireball 
effectively. The deceleration radius is Rd = l/T 2 J 3 . If A > l/T^/ 3 (the thick shell case), we 
define the deceleration radius as the shock crossing radius Rd = / 3 / 4 A 1 / 4 . The deceleration 
lab time is given by i d ~ Rd/ c- 

If the Blandford-McKee (BM) blast wave scaling V oc R^ 3 ^ 2 (Blandford and McKee 
1977) was valid for the whole fireball evolution (the thin solid line in figure 1), integrating 
eq. (1), we obtain R oc t 1//4 . Since the spectral characteristics of forward shock synchrotron 
emission are given by products of t, R(t), T(t) and constant parameters (Sari, Piran & 
Narayan 1998), the light curve should be described by a power-law with the reference time 
to = 0. However, the BM blast wave scaling is applicable only after the fireball is decelerated 
(e.g. Kobayashi, Piran & Sari 1999). At earlier times the fireball (the thick solid line) is in the 
coasting phase, and it is slower than evaluated from the BM blast wave scaling. The delay of 
photons from the shock front at the deceleration lab time is larger by 5t = (Sf — Sb) / c than 
in the case that the BM blast wave scaling is applicable to the whole evolution (see figure 
1) where Sb and Sp are the separations at the deceleration lab time between the gamma- 
ray front (dashed line) and the reference BM blast wave (the thin solid line), and between 
the gamma-ray front and the fireball forward shock (the thick solid line), respectively. The 
reference time for the afterglow modeling is given by 



In the thick shell case, a forward shock keeps being energized for a longer time, and the 
deceleration phase starts at a later time. It still overestimates the shift of the reference time 
if to is set at the end of the energy injection or equivalently at the peak time of afterglow. In 
figure 2 the dot indicates the point at which the energy injection stops and the fireball turns 
into a BM solution. Measuring times from the peak of afterglow corresponds to defining a 
reference null geodesic (the dashed line) as it goes through the dot. Imagine the reference BM 



1 In the thin shell case, the deceleration observed time is given roughly by the critical width ~ A c /c = 
Rd/cTl, and it is longer than the duration of the prompt emission A/c (Sari 1997). The deceleration time 
A c /c approaches A/c if we take a larger To, and around ro ~ T c = (l/A) 3 ^ 8 (or equivalently A ~ A c ) a 
reverse shock becomes relativistic during the propagation. The Lorentz factor of the shocked material at the 
crossing time becomes independent from the initial value, and it is given by r c . The deceleration observed 
time is about Rd/cT^. = A/c. 
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blast wave (the thin solid line) associated with the new reference null geodesic 2 . The decay 
of the emission from this reference blast wave should be characterized by a power-law using 
the new reference time (measuring times from the peak). The evolution of a real fireball (the 
thick solid line) is also described by a BM blast wave after the deceleration. However, it is 
a different solution (a different curve on the diagram). Applying the new reference time to 
the afterglow modeling leads to a wrong estimate of the decay index especially in the early 
phase. The decay index should become shallower than the real value. 

Using eq.(l) and the evolution of a BM blast wave T ~ (i?//)~ 3 / 2 , we obtain Sb = 
Rd/8T d where Yd is the Lorentz factor at the deceleration radius Rd- Since the Lorentz 
factor of a fireball is constant r = T at R < Rd, we get Sp = 4Ss for the thin shell. In the 
thick shell case, a reverse shock becomes relativistic before it crosses the fireball shell, and it 
begins to reduce considerably the Lorentz factor of the shell's matter which it crosses. The 
Lorentz factor of the forward shock is also significantly reduced as T ~ (I 3 / AR 2 ) 1 / 4 during 
the reverse shock crossing (Sari & Piran 1995; Kobayashi, Piran & Sari 1999), and we get 
Sf = 2Sb for the thick shell. The reference time is given by 



where td = Sp/c is the deceleration observed time. 

We have considered a simple broken power-law model for the evolution of a fireball 
Lorentz factor, and we obtained eq. (3). In reality the fireball should decelerate gradually 
around the deceleration radius. An artificial early steepening could happen only when Sb *C 
Sp, and in such a case to should be set at td (the afterglow peak time) to get the correct 
a. However, in both of the thin shell and thick shell cases, we have found that Sb and Sp 
are comparable. Therefore, measuring times from the beginning of the prompt phase to — 
should not induce an overestimate of afterglow decay right after the peak. 

Eq. (1) gives the delay of photons emitted from a point on the shock front on the line of 
sight, while most photons suffer longer delays, since they are emitted from a shocked region 
of finite thickness behind the shock, and from positions off the line of sight. Although these 
effects could make both separations S B and S F larger by a factor of a few (Waxman 1997), 
they are still comparable and our arguments are valid. 



2 BM blast wave lines on the spacetime diagram are determined by two parameters, an explosion time 
(i.e. a null geodesic on the diagram) and the Sedov length. The latter is evaluated with the total fireball 
energy if there is energy injection before the deceleration time. 




3td/4 the thin shell case 
t d /2 the thick shell case 



(3) 
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3. Numerical Model 

We employ a spherical relativistic Lagrangian code based on the Godunov method with 
an exact Riemann solver to evaluate the hydrodynamic evolution of a relativistic fireball 
(Kobayashi et al. 1999; Kobayashi & Sari 2000). Using the Einstein summation convention 
the equations describing the motion of a relativistic fluid are given by the five conservation 
laws 

d{pu l )/dx l = 0, 8T ik /dx k = 0, (4) 

where u l is four velocity and T lk is the stress-energy tensor (i,k = 0, ...,3), which for a 
perfect fluid can be written as T lk = wu l u k — pg tk . Here, g lk is the metric tensor, p the fluid 
pressure, and w = e + p the heat function per unit volume. Shocked material is extremely 
hot and pressure is related to internal energy e and mass density p as p = (e-pc 2 )/3. In the 
Godunov scheme, conservative variables are considered as piecewise constants over the mesh 
cells at each time step and the time evolution is determined by the solution of the Riemann 
problem (shock tube) at the inter-cell boundaries (e.g. Marti & Miiller 1999 and references 
therein) . 

Because of the relativistic beaming effect, the radiation from a jet before the jet break 
can be described by a spherical model with an isotropic energy. The initial configuration 
for our simulation is a static uniform fireball surrounded by uniform cold ambient material 
(ISM). It is determined by four parameters: an isotropic energy E, a dimensionless entropy 
r , an initial radius Ro and the ISM mass density p\. E and p\ always appear as the 
ration E/pi in the hydrodynamics computation, the system is actually determined by three 
parameters, the initial radius R , the entropy T and the Sedov length I. For convenience, 
we set the initial lab time as Rq/c rather than zero in numerical calculations, and the 
observed time is determined by the delay of photons from an emitter with respect to the 
photons emitted from the initial fireball surface at R = Rq and t = Rq/c. Since the fireball 
immediately accelerates to a relativistic velocity, this is practically equivalent to the observed 
time used in figure 1. 

We carry out numerical calculations for the total number 3 of mesh cells N = 540. 
A third of the cells are for the fireball, while the other cells describe the ambient medium 
within Rmax = 10 18 cm. Although the initial configuration is ~ 55 cells per decade, the 
nature of Lagrangian method gives a much higher resolution for shocks which sweeps cells 
and compresses them by a factor of ~ T 2 . A reflection boundary and free boundary condition 



3 We have reevaluated the light curves Figs 5 and 7 for N — 5400. In the log-log space, the resulting 
light curves are almost identical. The differences are less than several percent. When we remove numerical 
oscillations, they are in a few percent agreement. 
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are imposed at the center R = and at R ma x, respectively. We consider two cases - the thin 
and thick shell cases. 



3.1. The Thin Shell Case 

We first consider the thin shell case: E = 10 53 ergs, T = 100, R = 3 x 10 n cm, 
and pi = 1 m p cm" 3 where m p is the proton mass. Initially, as the fireball expands into a 
surrounding medium, a narrow shell with a radial width A ~ R is formed. The Lorentz 
factor of the shell increases linearly with the radius during the free acceleration stage. Then, 
the fireball shell uses up all its internal energy, and it coasts with the Lorentz factor of 
r . The coasting ends once the ISM begins to influence the shell. After the deceleration 
radius Rd, the profile of the shocked ISM medium begins to approach the BM solution. The 
evolution of a fireball is fully discussed by Kobayashi et al. (1999). 

The interaction between the shell and the ISM is described by two shocks: a forward 
shock propagating into the ISM and a reverse shock propagating into the shell. Figure 3 
shows the propagation of the shocks. Initially, the unshocked fireball shell has all the energy 
of the system. As the shell expands, the reverse shock decelerates the ejecta while the 
forward shock accelerates the ISM. The energy is transfered from the unshocked shell to the 
ISM via the shocks, finally the shocked ISM carries essentially all the energy of the system. 
In the intermediate stage, around the deceleration time t d = (3E /32ttpiC 5 Tq) 1 ^ 3 ~ 195 sec, 
the shocked shell has comparable energy to the shocked ISM. The evolution of the energies 
in three regions, inside of the reverse shock (unshocked shell), between the reverse shock 
and the contact discontinuity (shocked shell), and between the contact discontinuity and the 
forward shock (shocked ISM) are shown in Figure 4. The observed times of photons from 
the reverse shock front, the contact discontinuity and the forward shock front are used to 
describe the evolution of the energies. 

Even after the deceleration time, the reverse shocked shell (the thin dashed-dotted line) 
carries a significant fraction of the system energy for a long time, 20% at t ~ lOt^ and 10% 
at t ~ 45td. In the thin shell case, the reverse shock is Newtonian or sub-relativistic in the 
frame of the unshocked shell (the deceleration by the reverse shock is not significant as we 
can see in the top panel of figure 3). It does not heat the shocked region well. The reverse 
shocked region is already cold at the shock crossing time (the deceleration time), and the 
shocked shell carries the energy mainly in the form of the kinetic energy E oc V . Assuming 
a power-law decay Y oc R~ 9 , we obtain E oc T oc t~#/( 1+2 9) ~ t~ - 4 for g = 2.2 (Kobayashi & 
Sari 2000). 
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This slow energy transfer should lead to the round-off of an afterglow peak. To evaluate 
the afterglow light curve, we consider here a simple case in which the energy of the magnetic 
field remains a constant fraction of the internal energy B 2 oc p. The electron random 
Lorentz factor evolves as 7 m oc p/p after the shock heating where p and p are the pressure 
and density of a fluid element. The typical synchrotron frequency in the observer frame is 
v m oc Y^ m B. Since a relativistic shock totally ionizes material which it crosses, the total 
number of electrons in a fluid element (mesh cell) is N e = 4irR 2 dRTp/m p where dR is a 
cell width in the lab frame. The spectral power at the typical frequency from a cell is given 
by F Um oc N e YB. Assuming a power-law distribution of the electron random Lorentz factor 
with index p, the observed flux is 



where we have assumed that the observational band is located between the typical frequency 
v m and the cooling frequency v c . In early afterglow phase, with the typical parameters, X-ray 
and optical band satisfy this condition 4 for the forward shock and reverse shock emission, 
respectively. We evaluate the flux and observed time of photons from each shocked fluid 
element, and superimpose the emission to construct the forward and reverse shock light 
curves. The results are shown in figure 5. Although at later times the light curves approach 
the expected power-law decays L oc t^ 1 (forward shock; the thick solid line) and L oc t~ 2 
(reverse shock; the thin solid line), the peaks are rounded as we expected. 

At the deceleration time (and after that), the light crossing time of the forward shocked 
region is comparable to the observed time A/c ~ R/cT 2 . To evaluate the effect of the 
thickness of the forward shocked region, we assume that all the shocked electrons emit 
photons at the forward shock front. The resulting light curve peaks at an earlier time 5 
t ~ tj/2, but the shape around the peak itself is similar. We conclude that the round-off of 



4 For the typical microscopic parameters e e = 10 -2 and eb = 1CP 3 , the numerical results shown in this 
paper actually satisfy this condition in the time ranges of the plots ( Figs 5 and 7 ) where e e and €b are the 
fractions of shock energy given to magnetic field and electrons at a shock, respectively. Choosing a larger es 
could produce an additional break at t ~ 10 3 — 10 4 s especially in X-ray light curves due to the passage of the 
forward-shock cooling frequency in X-ray band. However, the steepening is small Aa=l/4 and it happens 
well after the deceleration of the fireball. The cooling break does not affect our discussions about the onset 
of afterglows. 

5 A peak time is given by the sum of a deceleration time and a fraction of the light crossing time of a 
shocked region. If all the shock energy is assumed to be radiated at the shock front, the emission should 
peak earlier. The analytic deceleration time td is obtained based on a simple model, and it could have an 
error of a factor of a few. 




(5) 
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the afterglow peak is mainly due to the slow energy transfer from a fireball shell to a blast 
wave. 

The internal shock model requires a highly irregular outflow from the GRB central 
engine. Since the hydrodynamic interaction inside the flow smooths the velocity and pressure 
profiles, but not the density profile, fireball ejecta might have an irregular density profile at 
the deceleration time. Emission from the ejecta during a reverse shock crossing could reflect 
the light curve of the prompt emission produced by internal shocks (Nakar & Piran 2004). 
The peak time of the reverse shock emission also depends on the density profile. As a fireball 
expands, a reverse shock accelerates more and more electrons in the fireball while the Lorentz 
factor and pressure of the shocked region decreases. The balance determines the peak time. 
In the thin shell case, a reverse shock does not effectively decelerate the shell material that it 
crosses. This leads to large pressure gradient across the shocked region. The contribution of 
emission from the inner part of the shell becomes negligible, and the reverse shock emission 
peaks slightly before the shock crossing time (the thin solid line) . 



3.2. The Thick Shell Case 

We consider the thick shell case in which a reverse shock becomes relativistic during 
the propagation. The initial condition is E — 10 52 ergs, T = 10 3 , R = 3 x 10 11 cm, and 
p 1 = 10 m p cm" 3 . The fireball is decelerated around R d ~ 1.5 x 10 cm, which 

corresponds to the shock crossing time td — Ro/c — 10 sec. 

The energy transfer from the fireball to a blast wave is similar to that described in the 
thin shell case. At the shock crossing time, the forward-shocked ISM and reverse-shocked 
shell have comparable energies (the thick lines in figure 4). The main difference is that in the 
thick shell case the shocked shell carries the energy in the form of the internal energy, instead 
of the kinetic energy, because the reverse shock significantly decelerate the shell. Figure 6 
depicts the profiles of T, p and p. When the reverse shock crosses the shell with a width 
of ~ Rq, a rarefaction wave begins to propagate towards the contact discontinuity, and it 
quickly transfers the shell's internal energy to the shocked ISM in a timescale comparable 
to the deceleration time (the shock crossing time). The steep decay of the shocked shell 
energy right after the peak (the thick dashed-dotted line) is due to the rarefaction wave 
propagation. Since for our parameters the reverse shock is mildly relativistic in the frame 
of the unshocked shell f ~ r /(//-Ro) 3 ^ 8 ~ several (Sari & Piran 1995), the shocked shell 
becomes cold after the propagation, and the energy of the shocked shell begins to decay in 
the same rate as the thin shell does. 
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In the thick shell case, the energy of a shocked shell is swiftly relayed to a forward shock 
by a rarefaction wave. A broken power-law T oc i?" 1 / 2 (R < R^) and Y oc R~ 3 ^ 2 (R > Rj) 
describes the evolution of the forward shock well. The round-off of an afterglow peak is 
expected to be small. Using eq. (5) and numerical results, we evaluate the light curves of 
the forward shock and reverse shock emission, which are shown in figure 7. Initially the 
luminosity of the both shocked region slowly increases as ~ t 1 / 2 (Kobayashi 2000), and they 
peak around the shock crossing time. The forward shock light curve (the thick solid curve) 
is slightly steeper right after the peak compared to at late times. Although the decay is 
described well by a single power law if the reference time to is set at the middle between 
the GRB trigger and the afterglow peak as we discussed in section 2, measuring times from 
the GRB trigger also provides a reasonable estimate for the afterglow decay. The reverse 
shock emission (the thick dashed curve) drops sharply a ~ 5 during the rarefaction wave 
propagation. The high latitude (off-axis) emission could contribute to the early flux (Kumar 
& Panaitescu 2000). The numerical fireball shell has a long tail in the density profile, and 
the shocked tail contributes to the light curve at late times. The emission from the outer 
part containing 80% of the shell mass gives a steeper light curve (the thin solid curve). 



4. Inhomogeneous Fireball 

If a burst has a precursor, we might need a shift of the reference time to a later time. 
A simple example is that a precursor with a negligible energy is followed by the main com- 
ponent. The quiescent period could be due to a dormant period in the central engine or due 
to the property of the outflow (e.g. homogeneous outflow in which prompt gamma-rays are 
not produced by internal shocks). In the former case, t should be set at the explosion time 
of the main component. If the separation between the precursor and the main component 
is much larger than the deceleration time of the main component, measuring times from 
the precursor introduces an artificial steeping in the early afterglow. Assuming a precursor 
and measuring times from it, we have replotted figures 5 and 7. If the precursor is located 
at At ~ 2000sec ~ 10t^ or ~ 3500sec ~ 18t^ before the main component, the decay index 
of the early forward shock emission could be overestimated as a ~ 3 or ~ 5 in the thin 
shell case. In the thick shell case, At = 7t d and 12t^ lead to steep decay indexes of a ~ 3 
and ~ 5, respectively. Since afterglow light curves are usually plotted with to setting at the 
GRB trigger time, precursors, although they might be energetically small, need to be strong 
enough to trigger GRB detectors (e.g. BAT) in order to cause an artificial steep decay. 
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4.1. Case Studies 

X-ray flares were originally reported from BeppoSAX observations, GRB 011121 and 
GRB 011211 (Piro et al 2005). GRB 011121 was the second-brightest GRB observed by 
BeppoSAX in gamma-rays (after GRB 990123) and in X-rays (after GRB 010222). The 
fluence in the 2-700 keV range corresponds to an isotropic energy of 2.8 x 10 52 ergs at the 
redshift of the burst z = 0.36. The gamma-ray light curve shows a main peak starting at 
t ~ 5 sec and ending at t ~ 30 sec with minor substructures. An X-ray flare took place at 
t = 240 - 310 sec. The fluence of GRB 011211 in the 2-700 keV range gives an isotropic 
energy of 3.6 x 10 52 ergs at the redshift of the burst z = 2.14. The gamma-ray prompt 
emission has a long duration T ~ 400 sec, and an X-ray flare is detected from 600 to 700 
sec. In the source frame, the two flares occurred at a similar time ~ 200 sec and they have a 
width A/c ~ 30 — 50 sec (see Piro et al. 2005 and references therein for the basic parameters 
of these bursts). 

Since the decay part of the X-ray flares and the following shallower part are described 
with a single power-law when the time is measured starting from the flare peak, the flares 
were suggested as the beginning of the afterglow caused by a thick shell (Piro et al. 2005). 
The thick shell interpretation requires a long energy release from the central engine. The 
burst should be long and extend all the way to the flare peak. However, the observed prompt 
emission ended well before the onset of the flares. In order to suggest the thick shell case, 
one must assume that the efficiency of conversion into gamma-ray varies dramatically with 
time. Furthermore, the engine should release most fireball energy at the last moment (major 
reenergization at the flare), because measuring times from the GRB trigger does not lead 
to an overestimate of the decay index even in the thick shell case, as long as the shell is 
homogeneous. 

In principle, the low conversion efficiency could originate from (1) a small Lorentz factor 
of late ejecta which does not allow gamma-ray radiation or (2) a small dispersion of the 
Lorentz factor of late ejecta in the internal shock scenario. The critical Lorentz factor is 
given by T c = (3£/47rp lC 2 A 3 ) 1/8 ~ 310E^ 8 A^ /8 where Q n = Q/10" in the cgs unit 
and A12 ~ A/(c x 40sec). In the scenario (1), the low Lorentz factor condition T < 100 
requires that the last energetic ejecta with an energy E > 10 53 ergs and a width A/c ~ 40 sec 
should be in the thin shell regime. The deceleration observed time t d > 420 E^ 3 p^ 3 T 2 8 ^ 3 
sec becomes much larger than the width of the flares ~ 40 sec, and one finds that the scenario 
(1) is inconsistent 6 . 



6 In the afterglow modeling of this thin shell case, measuring times from the GRB trigger does not cause 
an early steep decay, because the separation between the beginning of the prompt emission and the flare 



- 12 - 



To examine the scenario (2), we consider an inhomogeneous shell (two components) 
expanding into ISM. Corresponding to the major energy release from the central engine at 
the last moment, the inner edge of the shell with a width A^/c ~ 40 sec is assumed to 
have an energy E b larger than the energy E a ~ 10 53 ergs in the preceding outer part with a 
width A a /c ~ 200 sec. The Lorentz factors of the two components at the end of the internal 
shock phase are the same value of T. Both components should be in the thick shell regime, 
otherwise the deceleration time of the shell becomes larger than the flare occurrence time 
~ 200 sec. The deceleration radius of the shell should be larger than the radius R ~ A^ 4 /^ 4 
at which a reverse shock crosses the outer component where l a = l(E = E a ). The separation 
at the deceleration time between the gamma-ray front and the reference BM blast wave 
satisfies Sb > (A a /8)(E a / Ef,). The separation at that time between the gamma-ray front 
and the fireball is Sp ~ A a + A^/4 ~ A a where we have assumed the fireball evolution 
T oc R~ 1 / 2 before the deceleration. Note that a broke power law description of Y around 
the deceleration radius is a good approximation in the thick shell case. Then, we obtain 
Sf/Sb < 8(Et,/E a ). The precursor discussion at the beginning of this section corresponds 
to the two component model with E a = 0. We can show that the separation ratio at the 
deceleration time is S F /S B ~ 8(A a /A b ) in this case. Since replotting figure 7 we have found 
that A a /A fe > 10 to produce an early steep decay, in the scenario (2) E b should be at least 
ten times larger than E a . This requires a large energy budget for the central engine E > 10 54 
ergs. Even in the limit of E a /E b = 0, the small ratio A a /A fe ~ 200/40 = 5 does not lead to 
a very steep decay of a > 3 — 5. Therefore, we conclude that even an inhomogeneous fireball 
can not produce X-ray flares in early afterglow via external shock emission process. 



5. Conclusions 

We have studied the reference time to for the afterglow modeling. Although measuring 
times from the beginning of the prompt emission (GRB trigger) might cause a slight overes- 
timate of the early afterglow slope in the thick shell case. This choice of to gives a reasonable 
approximation, and it does not induce a very steep decay (a ~ 3 — 5) like the early steep 
decay or X-ray flares in the canonical Swift X-ray light curve. 

The leading model to explain the rapid decay and flares in early X-ray afterglow is the 
internal shock emission. A clear, testable prediction of this model is that the temporal decay 
index a of the tail part should be related to the spectral index (3 by an equation a = 2 + (3 



At ~ 200 sec is smaller than the deceleration time td > 420 sec. As we have shown at the beginning of this 
section, a significant artificial steepening happens only when the ratio At/td is larger than ~ 10. 
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(Kumar & Panaitescu 2000). When evaluating the emission decay L oc (t — t )~ a in the 
internal shock model, an important difference is that the GRB trigger time is no longer 
special. The reference time to should correspond to the onset of each particular spike in 
the prompt emission or in afterglow (Kobayashi et al. 1997; Zhang et al. 2006a; Nousek et 
al. 2006; Fan & Wei 2005). Every time when the central engine is re-stared to eject sub- 
shells, the reference time to should be re-set to the reactivation time of the engine. Although 
O'Brien et al. (2006) have found that a appears to be largely independent of (3 when the 
BAT trigger time is used as t , Liang et al (2006) have shown that the relation a = 2 + (3 
is more or less satisfied in most cases if to is set near the beginning of rising segment of 
the last pulse of the prompt emission or a corresponding X-ray flare, and if the underlying 
forward shock emission component is subtracted. Swift observations support the internal 
shock model. 

The self-consistent internal shock interpretation should be more favorable than the 
beginning-of-the-afterglow interpretation. The latter can not explain multiple X-ray flares 
in a single event. Such behavior is observed in many Swift bursts (Burrows et al. 2005; 
Falcone et al. 2006; Romano et al. 2006; O'Brien et al 2006). Usually X-ray flux already 
begins to decay before X-ray flares appear, and it suggests that the onset of afterglow is 
prior to the flares. If a large amount of energy is impulsively injected to a fireball during the 
deceleration, to might be re-set to the injection time. However, an afterglow baseline also 
should shift after a flare (energy injection). This clearly contradicts with observations in 
which after a flare peak, afteglow decays back to its pre-flare flux level. We can not explain 
X-ray flares by the shift of t associated with large energy injections. 

Swift discovered that a large fraction of X-ray afterglows have a slow decay phase, and 
it is suggested that energy injection into a blast wave takes place several hundred seconds 
after the burst. This implies that right after the burst the kinetic energy of a blast wave 
is very low and in turn the efficiency of internal shock process is extremely high (Zhang et 
al. 2006a; Nousek et al. 2006; Ioka et al. 2005; Granot, Konigl & Piran 2006; Zhang et 
al. 2006b; however see also Fan & Piran 2006). The round-off forward shock peak in the thin 
shell case might be a good candidate for the shallow decay phase. However, if we interpret 
the observed shallow decay as the round-shape of the deceleration phase, the model light 
curve is shallower than some of the observed ones. It may be because the observed curve is 
the combination of this round phase and the rapid decay from the GRB tail emission (high 
latitude). Equalizing the deceleration time and the shallow phase timescale we obtain the 
initial Lorentz factor r ~ 110( 3 / 8 E^ 8 p^^t^ 3 / 8 where the time dilation effect is taken into 
account ( — (1 + z)/3. 

For a wind environment, we can discuss the to issue in a very similar way. The Lorentz 
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factor of a shocked shell is constant in both of the thin and thick shell cases during the 
shock crossing;, while the BM blast wave decelerates as T oc R- 1 / 2 (Kobayashi & Zhang 
2003). The separations are comparable Sp ~ 2>Sb at the deceleration time. Measuring 
times from the beginning of the prompt phase should be a good approximation for events 
in a wind environment also. Since most bursts in a wind environment fall in the thick shell 
case (Kobayashi, Meszaros & Zhang 2004), and since the deceleration time is comparable to 
GRB duration in the thick shell case, another process (e.g. refreshed shocks) rather than 
the afterglow peak is necessary to explain the shallow decay phase. 

We thank Luigi Piro for useful discussion. This work is supported by NASA NNG05GB67G 
and NNG06GH62G. 
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Fig. 1. — Spacetime diagram: If the evolution of a fireball satisfies the BM blast wave scaling 
T oc i?~ 3 / 2 from the beginning, the thin solid curve gives the trajectory. However, the fireball 
initially coasts with a finite Lorentz factor. The evolution is described by the BM solution 
only after it begins to be decelerated. At the deceleration lab time the fireball (thick solid) 
gets behind the reference BM blast wave (thin solid), and the trajectories approaches each 
other at a later time. The observed time is given by the delay of photons from an emitter at 
a lab time with respect to the gamma-ray front (null geodesic; dashed line). Null geodesies 
are trajectories of photons. Physical trajectories should be steeper or parallel to this minimal 
slope line. A trajectory of an object at rest is vertical. 
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Fig. 2. — wrong choice of t : Fireball Evolution (thick solid curve), a shifted null geodesic 
(dashed line) and its reference BM blast wave (thin solid). 
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Fig. 3. — Thin Shell Case: Profiles of T, p and p at different lab times i = O.Q2R d /c, 
0.8QR d /c, l.lAR d /c, lA5R d /c and 1.74^/0. The x— axis is the distance from the contact 
discontinuity (CD). Photons emitted from the discontinuity at the lab times are observed 
at t = 100,200,400,840 and 1650 s, respectively. The gamma-ray front is located at (R — 
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Fig. 4. — Energy transfer from the fireball to the shocked regions, the thin shell case (thin 
lines) and the thick shell case (thick lines). The sum of the kinetic energy and the thermal 
energies in the unshocked shell (dashed lines), in the reverse shocked shell (dashed dotted 
lines) and in the forward shocked ISM (solid lines). All the energies are normalized by the 
explosion energies. The observed time is normalized by the deceleration time id- 
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Fig. 5.— Afterglow Light Curve: the Thin Shell Case. E = 10 53 ergs, T = 100, R = 
3 x 10 11 cm, px — 1 m p cm" 3 and p = 2.2. Forward shock emission in X-ray band (thick 
curve) and reverse shock emission in optical band (thin curve). Power-law fits (dashed lines): 
forward shock a — 1.0 and reverse shock a = 2.1. The light curves are normalized at the 
peaks. 
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Fig. 6. — Thick Shell Case: Profiles of T, p and p at different lab times t = 0.75R d /c, 
lA0R d /c, lA5R d /c, 1.92R d /c and 2AAR d /c. The x— axis is the distance from the contact 
discontinuity (CD). Photons emitted from the discontinuity at the lab times are observed at 
t — 5, 18, 19, 50 and 100 s, respectively. The gamma-ray front is located at (R — Rcd) = ct. 
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Fig. 7.— Afterglow Light Curve: the Thick Shell Case. E = 10 52 ergs, T = 1000, R = 
3 x 10 11 cm, pi = 10 m p cm -3 and p = 2.2. Forward shock emission in X-ray band (thick 
solid curve) and reverse shock emission in optical band (thick dashed curve). Reverse shock 
emission from the outer part of the shell which corresponds to 80% of the shell mass (thin 
solid curve). Forward shock emission decays slightly faster at the beginning as t~ 13 and later 
it decays as t~ 10 . reverse shock decay indexes a = 1.5 and a = 1.7. The deceleration time 
td — Rq/c — 10 sec (vertical dotted line). The light curves are normalized at the peaks. 



